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ABSTRACT 

Context. Interpretation of light curves of many types of binary stars requires the inclusion of the (cor)relation between surface 
brightness and local effective gravity. Until recently, this correlation has always been modeled by a power law relating the flux or the 
effective temperature and the effective gravity, namely Teff 

Aims. We look for a simple model that can describe the variations of the flux at the surface of stars belonging to a binary system. 
Methods. This model assumes that the energy flux is a divergence-free vector anti-parallel to the effective gravity. The effective 
gravity is computed from the Roche model. 

Results. After explaining in a simple manner the old result of Lucy (1967, Zeit. fiir Astrophys. 65,89), which says that /? ~ 0.08 
for solar type stars, we first argue that one-dimensional models should no longer be used to evaluate gravity darkening laws. We 
compute the correlation between log T^ff and log geff using a new approach that is valid for synchronous, weakly magnetized, weakly 
irradiated binaries. We show that this correlation is approximately linear, validating the use of a power law relation between effective 
temperature and effective gravity as a first approximation. We further show that the exponent fj of this power law is a slowly varying 
function, which we tabulate, of the mass ratio of the binary star and the Roche lobe filling factor of the stars of the system. The 
exponent p remains mostly in the interval [0.20, 0.25] if extreme mass ratios are eliminated. 

Conclusions. For binary stars that are synchronous, weakly magnetized and weakly irradiated, the gravity darkening exponent is well 
constrained and may be removed from the free parameters of the models. 
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1. Introduction 

Gravity darkening refers to the variations of the energy flux at 
the surface of a star that result from the variations of the local 
effective gravity. This eff'ect is especially visible on rapidly ro- 
tating stars leading to a dark equator and bright poles. It is now a 
key effect to interpret the brightness distribution of rapidly rotat- 
ing stars when observed by interferometry and to derive the right 
posi tioning of the rotation axis with respect to the line of sight 
(e.g. lDomiciano de Souza et al.ll2005l : iMonnier et al ] l2007l) . 

However, gravity darkening also affects the tidally distorted 
stars in multiple systems. It has long been known that the inter- 
pretation of the light curves of eclipsing binaries, in particular 
semi-detached or contact ones, requires the inc lusion of grav- 
ity darkening (e.g. iDiurasevic et al.ll2003l l2006h . The so-called 
ellipsoidal variables are understood as binary stars whose light 
variations come fro m the rotation of their tidally distorted figure 
dWilson et al.ll2009h . 

The theoretical modelling of gravity darkening started long 
ago with the work of von Zeipel ( 1924) who discovered that in 
barotropic stars the energy flux is proportional to the local effec- 
tive gravity leading to the now celebrated von Zeipel law which 
relates effective temperature and effective gravity by Teff oc g^^^. 
However, this law also assumes a radiative envelope for the star. 
Thus, motivated by the importance of grav ity da r kening in con- 
tact binaries that harbour low mass stars, ILucvI (119671) studied 
how to generalize this law to stars with a convective envelope. 
Still considering a power-law relation, Tefj oc and calibrating 
it on one-dimensional models, Lucy inferred that yS ^ 0.08 for 
solar- type stars. 



However, the importance of gravity darkening in stellar 
physics motivated further investigations of the theo retical prob- 
lem to give better predictions of the exponent yS (e.g. Clared 19981 
2003, 2012). But these works, as Lucy's one, are based on ID - 
models of stellar structure and lead to gravity d arkening lavv s 
that depend on the det ails of the atmosphere (e.g.lClarejl2012h . 
Moreover, as noted in lEspinosa Lara & Rieutordl (1201 ih . these 
laws are valid for small deviations from spherical symmetry. It is 
therefore not surprising that some observations do not fit because 
observed gravity darkening is generally detected on objects that 
show a strong effect and are far from spherical. 



In lEspinosa Lara & Rieutordl ( 1201 lb . hereafter referred to as 
paper I, we derived a model of gravity darkening for rotating 
stars that is not restricted to small deviations from sphericity. 
Moreover, it has no free parameter and perfectly matches re- 
cent observational results on rapidly rotating stars like a Leo 
or a Aql. This model rests on very few assumptions, essentially 
that the flux is anti-parallel to effective gravity, and should be 
very useful for studies of tidally distorted stars. 



In this paper we generalize the work of paper I to tidally 
distorted stars in order to provide a robust estimate of the corre- 
lation between flux and local surface gravity. This is achieved in 
sections 3 and 4. Before that, we feel that we ought to discuss 
the general problem of gravity darkening in spheroidal stars and 
the previous attempts to solve it (sect. 2). 
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2. The problem of gravity darkening laws in 
spheroidal stars 

2.1. Introduction 

Before getting into the details of modelling gravity darkening 
in tidally distorted stars, we reconsider the general problem of 
deriving a gravity darkening law, namely a law that would relate 
the local energy flux and the local effective gravity. 

At the surface of a star the energy flux and eff'ective gravity 
are two vectors that depend on the colatitude 6 and longitude ip: 

F = F(0,ip) and geff = g,d0,v) 

The existence of a universal relation between these two quanti- 
ties is not obvious because they are not locally related by any 
physical process (they are not local thermodynamic al state vari- 
ables!). However, if we consider the star to be barotropic, that 
is that the stellar fluid obeys an equation of state where pressure 
solely depends on density, namely P = P(p), then all thermo- 
dynamical variables only depend on the total potential (gravi- 
tational plus centrifugal) This implies that the radiative flux 
F - -Xr^T is proportional to g-eff = -VO because: 

F = -Xri^Wrm = -;t',(<l))r'(0)V<D ^XrmT'mg^s (1) 

While the surface of a star is generally defined by some value of 
the opacity, the surface of a baroropic star can be equivalently 
defined by an equipotential (O - (tj). The foregoing relation 
shows that the surface variations of the flux f just follow those 
of the effective gravity, which is the famous von Zeipel law. 

The assumption of a barotropic fluid is strong because it is 
enforced locally. However, the actual thermodynamic structure 
of a star never deviates strongly from barotropicity : in convective 
region matter is strongly mixed, which homogenizes entropy, 
while in radiative regio ns the recent two-dimensional mo dels of 
rapidly rotating stars of 'Es pinosa Lara & RieutordI (1201 1 ) show 
that the relative inclinations of isotherms, isobars or equipoten- 
tial is less than half a degree. 

In this perspective, the result of iLucvl (Il967h . who predicts 
jS ^ 0.08, which is a weak dependence of the flux with local 
gravity, is quite strange. Moreover, it should apply to stars with 
a convective envelope, which is an envelope close to isentropy. 



2.2. Flux-gravity law from spherically symmetric envelopes 

In order to clarify this situation, it is worth trying to understand 
the origin of Lucy's result. 

Let us first assume that in the range of temperature and den- 
sity of the envelope, the opacity of the matter can be approxi- 
mated by a power law of the form; 

K = KoffT-' (2) 

It was noticed bv lChandrasekhaii ( ll939h that a spherically sym- 
metric radiative envelope where both self-gravity and nuclear 
heat generation can be neglected behaves like a polytropic gas 
of index 



Surface pressure 



n = 



s + 3 

;7TT 



(3) 



At the surface of the star, one dimensional models usually as- 
sume that the star radiates like a black body while it is in hy- 
drostatic equilibrium. These hypothesis are condensed into the 
boundary condition 



3 K 



Gra vity 



High 
pressure 



high 

temperaiure 



Low 
pressure 



Fig. 1. Schematic representation of the generation of a baroclinic 
flow: the two columns represent two columns of fluid staying 
at two different latitudes; the difference in the temperature pro- 
file and the assumption of hydrostatic equilibrium in the column 
leads to two different pressure at their base: the fluid flows then 
from the high pressure region to the low-pressure one. This cir- 
culation transports heat in the latitudinal direction. In the rotat- 
ing background of a star, the Coriolis acceleration transforms 
this flow into a differential rotation (the so-called thermal wind 
in the Earth atmosphere). 



whi ch is largely used in stellar evolution code s (e.g. 
Kip penhahn & WeigertllT998t iHansen & Kawaledll994 . In this 
expression g is the surface gravity, P the surface pressure and k 
an averaged opacity. We note that this condition is used for both 
types of envelopes, convective or radiative. This may be surpris- 
ing in the case of a convective envelope but we recall that in the 
uppermost layers of a convective zone, convection is very ineffi- 
cient and flux is almost only radiative. 
Since the atmosphere is polytropic 



Pec r" 



and 



Thus we find: 

g oc 

leading to the gravity darkening exponent: 



1 



n{^i +1) + 1 - s 



(5) 



(4) 



We immediately observe that the expression of the polytropic 
index (|3]l implies that for a radiative envelope 

1 

showing that this exponent is actually independent of the expo- 
nents of the opacity and gives back the von Zeipel law, as ex- 
pected. 

We now turn to the convective envelope. In this case the fluid 
is assumed almost isentropic and therefore the relation between 
thermodynamic variables P, p, r is enforced by the « = 3/2 poly- 
tropic index, which therefore no longer depends on the opacity. 
However, the condition (HJl is still enforced and the opacity may 
still be described (locally) by power laws. Thus expression (|5]l 
can be used with a polytropic exponent n - 312 and the local 
exponent of the opacity laws. At the surface of the Sun, the fit- 
ting formula proposed by J. Christensen-Dalsgaard for simple 
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solar models (e.g. IChristensen-Dalsgaard & Reiteilll995h . sug- 
gests = 0.407895 and s = -9.28289 thus leading to 

/3 ^ 0.0807 



which is just the value found bv lLucvl (Il967h for solar type stars. 

This result therefore explains the low values that have been 
derived from one-dimensional spherical models of solar type 
stars and shows why the exponent is sensitive to the details of 
the opacity in the envelope. 

2.3. The case of non-spherical stars 

Until recently, the lack of multi-dimensionnal models able to 
describe stars distorted by rotation or tidal forces, lead to grav- 
ity da rkening laws derived from one-dimensional models (Claret 
Il998i 120121) . This expedient is built on the assumption that (|4]l 
is still true in these spheroidal stars, which is dubious since (|4]l 
rests on the hydrostatic equilibrium, namely 

dP^g 
dr K 

In rotating or tidally distorted stars, the mechanical balance near 
the surface is rather 

V Vv = --VP - VO 

P 

where v is the fluid velocity in some inertial frame. The new ve- 
locity terms disappear at the pole of the star but are important 
in the equatorial region. Thus, the mechanical balance now in- 
cludes the v elocity fields that are generated by the barochni c 
torque (e.g. RieutordI l2006t lEspinosa Lar a & Rieutord' '2007). 
The juxtaposition of ID-solutions is not mechanically possible 
because it generates a latitudinal pressure desequilibrium (Fig.[T] 
sketches out the mechanism leading to baroclinic flows). The 
derivation of a gravity darkening law should therefore take into 
account the bi- or tri-dimensionality of the models. 



3. The model 

3.1. introduction 

In order to circumvent the use of one-dimensional models, we 
propose to consider the constraints met by the flux in rotating or 
tidally distorted stars. First, we know that in the envelope of a 
star in a steady state 



V-f = 



(6) 



i.e. the energy flux F is conserved. Second, we may assume 
that the flux is almost anti-parallel to the local gravity. We have 
checked on 2D models that this assumption is likely very close 
to reality in isolated rotating stars, even in the most extreme case. 
In paper I, we have shown that the angle between the two vectors 
does not exceed half a degree, which is negligible. 

In tidally distorted stars, synchronously corotating with the 
orbital motion, this may well be the case because rotation is usu- 
ally weaker than in a star rotating near criticality. However, irra- 
diance by the companion may break this assumption. Obviously, 
the additional flux from the outside may have a dramatic influ- 
ence if strong enough. Thus, the following theory applies to bi- 
nary stars when the companion is weak enough not to perturb 
the thermal balance of the primary; other constraints possibly 



affecting low-mass stars are discussed in the conclusion section. 
To further simplify the problem we will also assume that the 
mass distribution of the primary is like a point mass. This is the 
usual Roche model, which is not a strong approximation since 
we are focussing on the gravity field in the envelopes. 

Provided these hypothesis, we can derive the dependence of 
the flux and the effective gravity with latitude and longitude. As 
we shall see these quantities are correlated (but not related) as 
expected from the previous considerations. To make this cor- 
relation easily usable, we will approximate it by a power law, 
namely Peff gffj, but this has no physical significance except 
of showing the importance of the discrepancy with von Zeipel 
law. 



3.2. Equations 

The main hypothesis is the same described in paper I. We con- 
sider that the energy is transported predominantly in the vertical 
direction indicated by the local effective gravity geif , that is, the 
flux vector F is always anti-parallel to the effective gravity 



(7) 



where / is an unknown positive function of the position to be de- 
termined and r, and tp are the spherical coordinates. Combining 
this relation with flux conservation (|6|l, we get 



(8) 



As in paper I, we will consider for simplicity that the star is 
centrally condensed and adopt the Roche model for the effective 
gravity. Then, all the energy is generated at the centre of the star, 
which implies that 



limf : 



/() 



ArGM 



(9) 



If we denote by n the unit vector in the vertical direction, such 
fliat getf - -ge«n, then we can write 



«■ Vlog/ 



(10) 



Eq. ([Tol l gives the variation of / in the local vertical direction. 
Let us build a trajectory C(6o, (fiu) that starts at the centre of the 
star with an initial direction indicated by (Oq, ipo) and is tangent to 
n at every point. C(0o, ^o) is therefore a field line of the effective 
gravity field. The value of / at a point r along the curve can be 
calculated as a line integral 

m = /„ exp I I dl] for r e C(0„, ^o) ■ (1 D 



/ 

JC{9o,<, 



According to the Roche model, the effective potential of a 
synchronized binary system with circular orbits is given by 



GM, 



GM2 



1 



Vfl^ + - 2ar cos 



- -Q^r^(sin^ 6 sin^ w + cos^ 6) + a 
2 ^ M1+M2 



M2 



(12) 



Q r cos , 



where a is the distance between the centres of the two stars and 
Ml and M2 are the respective masses. The coordinate system 
whose origin is at the centre of the primary star, is oriented so 
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Fig. 2. Schematic representation of the primary star with filHng 
factor 0.8. The position of the Lagrange point Liis indicated. 



that the z-axis corresponds to the axis connecting the stars cen- 
tres. The plane Oyz is the orbital plane and the ;c-axis is normal 
to this plane. This implies that the spherical angles are not the 
usual colatitude and longitude of the star. The orbital angular 
velocity is 



Q 



G(Mi + M2) 



We normalise the expression of the potential by GMi/a, so 
that 



1 



■-(1 + q)r(sm^ Osin^ ip + cos^ Q) + qrcosO , 



(14) 



where ^ - is the normalised potential and T = r/a, the 

normahsed radius. We have also introduced the mass ratio 



M2 
M, 



(15) 



Note that the normalised angular velocity is oj - ^fV^Tq. The 
effective gravity is ggff - -V0 and satisfies 



V-g,ff = 2(1+^) 



3.3. Solution of the model 



(16) 



We limit ourselves to the case of detached or semi-detached bi- 
naries, focusing on one of the stars that we call the primary. 
Unfortunately, eq. ( fTTT i has resisted all our efforts of finding a 
solution in a closed form forcing us to solve it by numerical in- 
tegration in the following way. 

Let ^ be the z-coordinate of the point at which an equipo- 
tential surface 4> =const. intersects the axis connecting the two 
stars, normalised by a the position of the secondary. Its maxi- 
mum value ^max corresponds to the surface of the primary. If ^1 
is the distance from the centre of the primary to the Lagrange 
point Li we define the filling factor as p = ^max/^i (see Fig.|2]). 
A star with p = 1 fills its Roche lobe. This definition is useful 



because, apart from normalising constants, the solution depends 
only on the mass ratio q and the filling factor p. 

As ^ marks a single equipotential, the potential can be ex- 
pressed as a function of ^ only, (p - (p(^) and 



(17) 



where goiO = Wgefiix = 0,y = 0,z = ^)||. 

We can use ^ as the parameter describing the field line 
C{Oo, (fio) that appears in (fTTl i. The differential equation of the 
field line is given by 



geff X dr = . 



(18) 



which expresses that the differential vector is parallel to the ef- 
fective gravity. Then, using ( fTTb and ( fTSl ), we can write 



dr 



^0 

— ^eff 



(19) 



8efl 



from which we can calculate the field line points r(f , Oq, (fo), for 
a given starting direction (Oq, (po). The line element is 



dl = lldrll = ^d^ , 

^eff 



(20) 



Combining (fTTT i. ( fTSl l and (l20l l. we get the function / which 
relates the energy flux and the effective gravity. We find 



(13) /(f^. '^^) = /o exp 



Jo 



2(1 +q) 



[^eif(r(^,0O,^o))] 



(21) 



where 6^ and (ps are the coordinates of the point where the field 
line C intersects the surface of the primary. 

Equations ( fTOl i and ( fT9] ) form a system of four differential 
equations for solving / that is well-suited for numerical integra- 
tion. Given an initial direction (0o, (/jq) and using initial condition 

the system can be integrated from the centre to the surface. 
We thus get the value of / at the surface for each {Oq, (f^) and 
thus for each {6^, (ps)- There is a one-to-one correspondance be- 
tween these coordinates because these field lines do not cross 
and effective gravity is defined everywhere. 

The integration has been performed using the Runge-Kutta- 
Fehlberg method, often referred to as RK F45, which combi nes 
two RK methods of fourth and fifth order (iPress et al.ll2002h . It 
provides an estimate of the truncation error that can be used to 
adapt the step size. 

4. Results 

Figure |3] shows the results for several combinations of the mass 
ratio q and filling factor p. The normalised values of effective 
temperature and gravity that appear in the figure are defined as 



T'pfF — 



/ 47ro-fl^ 



2\l/4 



!eff 



and 



?eff = 



GM, 



?eff 



(22) 



In the left, we can see the distribution of T^ff at the surface in one 
hemisphere of the star, projected in the plane XZ (remember the 
definition of the coordinate system shown in Fig.|2]). A logarith- 
mic plot of Teff versus geff for each case can be seen on the right. 
Since the problem has no azimuthal symmetry, the same value 
of geff appears at various places on the surface, but associated 
with different values of T^ff, thus illustrating the absence of a 
one-to-one relation between T^ff and g^ff. However, as shown by 
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Fig. 3. Left: Distribution of Teff at the surface in one hemisphere of the primary, projected in the plane XZ for several combinations 
of the mass ratio q - M^IMx and the filling factor p. Right: Logarithmic plot of Tetj versus geft for the same cases. Solid lines are 
linear fits to the data and the corresponding value of the gravity darkening exponent /3 is indicated in the plots. For comparison, von 
Zeipel's law (j3 = 0.25) is shown by a dashed line. In the second row an inset shows the small dispersion of the points around the 
fitting power law. The last case (q - 0) corresponds to an isolated star rotating with angular velocity Q = 0.854ilj:. Normalised 
values of effective temperature and gravity are defined in Eq.l22l 



the plots, the dispersion is very small (see the inset in Fig.[3]l. It 
emphasises a good correlation between the two foregoing quan- 
tities. Moreover, the log-log plot shows that this coiTelation is 
well approximated by a power law. The actual value of the expo- 
nent /3 is always slightly less than 1/4. The worst case is the last 



one where q - O.ln fact, this case coiTesponds to the limit of an 
isolated star (M2 = 0) rotating at angular velocity 

D. = p^'^ilk , (23) 
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Table 1. Fitted values of the gravity darkening exponent /? in the correlation between T^g and geft as a function of the mass ratio 
q - Mil Ml and the filling factor p. As the relation does not follow exactly a power law, the maximum relative error in T^s calculated 
using the fit compared with the exact value is included in parentheses. 



Mass ratio {q) 


p=0.3 


p=0.5 


p=0.7 


Filling factor (p) 
p=0.8 


p=0.9 


p=0.95 


p=0.99 


q=Ql.O 


0.2456 


0.2308 


0.2034 


0.1854 


0.1641 


0.1509 


0.1365 




(0.00004%) 


(0.0038%) 


(0.078%) 


(0.27%) 


(0.96%) 


(2.08%) 


(6.08%) 


9=0.001 


0.2464 


0.2342 


0.2110 


0.1958 


0.1789 


0.1708 


0.1666 




(0.00002%) 


(0.0021%) 


(0.043%) 


(0.15%) 


(0.48%) 


(0.92%) 


(1.55%) 


9=0.01 


0.2472 


0.2374 


0.2188 


0.2068 


0.1948 


0.1901 


0.1880 




(0.00002%) 


(0.0013%) 


(0.025%) 


(0.090%) 


(0.31%) 


(0.56%) 


(1.31%) 


9=0.1 


0.2482 


0.2420 


0.2303 


0.2232 


0.2168 


0.2143 


0.2131 




(0.0U0U47C) 


(U.UOlzyo) 


(0.0147c) 


/n r\/i nci \ 


(0. love) 


(U.JOTfo) 


(U.DOTc) 


9=1.0 


0.2489 


0.2452 


0.2385 


0.2346 


0.2308 


0.2291 


0.2280 




(0.00004%) 


(0.0013%) 


(0.011%) 


(0.050%) 


(0.23%) 


(0.52%) 


(1.44%) 


9=3.0 


0.2490 


0.2455 


0.2394 


0.2359 


0.2324 


0.2308 


0.2297 




(0.00004%) 


(0.0015%) 


(0.020%) 


(0.092%) 


(0.37%) 


(0.80%) 


(2.18%) 


9=10.0 


0.2489 


0.2453 


0.2391 


0.2355 


0.2320 


0.2304 


0.2292 




(0.00005%) 


(0.0019%) 


(0.036%) 


(0.15%) 


(0.55%) 


(1.16%) 


(3.09%) 


9=100.0 


0.2487 


0.2447 


0.2377 


0.2337 


0.2298 


0.2279 


0.2267 




(0.00005%) 


(0.0029%) 


(0.061%) 


(0.23%) 


(0.82%) 


(1.75%) 


(4.57%) 



where is the Keplerian angular velocity at the equator 
Qi = for which we recover the results obtained in pa- 

per I. 

In Table[T|we summarize the fitted values of the gravity dark- 
ening exponent as a function of q and p. As pointed out above, 
the relation between T^g an ^eft is not exactly a power law so, 
in order to measure the quality of the fit, we also included, in 
parentheses, the maximum relative error in Tgff calculated using 
the fit compared to the exact value. Again, for intermediate q, the 
departures from von Zeipel's value are moderate, with a value of 
yS above ~ 0.21 and a correlation between logTeff and loggeff 
that remains more or less linear. 

The largest deviation is observed for small q, for which the 
effects of nonlinearity are also the strongest. The same behavior 
can be observed in Fig. |4l where we plot the fitted values of ji 
versus the filling factor p for different values of q. As expected, 
the asymptotic limit of p = 0, i.e. spherical stars, corresponds to 
von Zeipel's value fi - 1/4. For a given p, the maximum value 
of j6 is attained for ^ ~ 3. 

5. Conclusions 

We have explained in a simple way the old result of lLucyl ( ll967h . 
who derived a gravity darkening law with an exponent fi ~ 0.08 
from ID models of solar-type stars. We trace this value of /3 as 
a consequence of the boundary conditions used by these stellar 
models, namely that surface pressure varies like g/k and that 
convective envelopes are close to a n = 3/2 poly trope. 

Noting that 2D or 3D models of rotating or tidally distorted 
stars cannot be in hydrostatic equilibrium, we argue that the use 
of gravity darkening exponents derived from ID models should 
be avoided. Following paper I, we propose a simple model of the 
gravity darkening effect that rests on hypothesis generally valid 
in the envelope of stars: 

1. flux is conserved (no heat sources); 

2. flux is antiparallel to local effective gravity; 

3. gravity is close to that of the Roche model. 
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Fig. 4. Fitted values of the gravity darkening exponent as a func- 
tion of the filling factor p for different values of the mass ratio 
q = M2/M1. Von Zeipel's value /3 - 0.25 is only attained in the 
limit p = (spherical stars). 

These three hypothesis lead to a correlation between effec- 
tive temperature and effective surface gravity which can be de- 
scribed by a power law with a gravity darkening exponent /3 that 
weakly depends on mass ratio and Roche lobe filling factor If we 
discard extreme mass ratios, /? remains in the interval [0.2, 0.25]. 
It is also independent of the nature (convective or radiative) of 
the envelope. 

This result is interesting to estimate gravity darkening in 
bright (hot) components of binary stars. It shows that the as- 
sumption of the von Zeipel value /3 - 1 /4 is certainly too restric- 
tive. 

However, a blind use of this model in binary stars is not rec- 
ommended either Indeed, this model does not include: 

a. irradiation from the secondary. 
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b. Coriolis effects on convection, 

c. magnetic fields, 

d. possible non-synchronism. 

The first two effects may indeed invalidate the assumption of 
the antiparallelism between flux and gravity. Magnetic field may 
generate spotted regions thus destroying the correlation between 
effective gravity and flux. Finally, non-synchronism, which may 
come from differential rotation, alters the assumed shape of the 
star 

The first three effects may especially affect low-mass com- 
panions of binary stars. For these stars, it is likely that the cor- 
relation between surface effectiv e gravity and effective temper- 
ature is perturbed. The results of jPiurasevic et all (I2003L l2006l) 
on semi-detached binaries, which give small gravity darkening 
exponents (namely p ~ 0.1) for late-type companions, may be 
viewed as the signature of a weak correlation between T^s and 
geff for these stars. 

In conclusion, the gravity darkening exponent should be in- 
terpreted as the parametrization of a correlation between T^s and 
^etf • For binary star components, where conditions (a-b-c-d) do 
not apply, this exponent should be fixed to a value derived from 
Table 1 according to the mass-ratio of the system and the fill- 
ing factors of the stars. If one ore more effect of the list exists, 
clearly a more elaborated model is needed. 
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